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Abstract. We present the generalization of the Sedov- Taylor self-similar strong 
spherical shock solution for the case of a central energy source varying in time, 
E — At k , where A and k are constants. The known Sedov- Taylor solution corre- 
sponds to a particular adiabatic case of k — or instant shock with an instant 
energy source of the shock, E — A. The self-similar hydrodynamic flow in the 
nonadiabatic k ^ case exists only under the appropriate local entropy (energy) 
input which must be supported by some radiative mechanism from the central 
engine. The specific case of k = 1 corresponds to a permanent energy injection 
into the shock, or injection shock with a central source of constant luminosity, 
L = A, E = At. The generalized self-similar shock solution may be applied to 
astrophysical objects in which the duration of central source activity is longer 
than the shock expansion time, e. g. the early phase of SN explosions, strong wind 
from stars and young pulsars, non-steady spherical outflow from black holes and 
collapsing dense stellar clusters with numerous neutron star collisions. 

Key words, hydrodynamics - shock waves - stars: winds, outflows - supernovae: 
general - ISM: bubbles 



1. Introduction 



The well known spherical shock solution (Sedov 1946, von Neumann 1947, Taylor 195C 



Stanyukovich 1969] ) describes the self-similar expansion of a strong spherical shock gen- 



erated by instant deposition of energy E — const by the central source in a homogeneous 
gas medium with density p\ = const. This instant shock solution is commonly used in 
numerous astrophysics applications e.g. for the modeling of SN explosions and evolution 
of young SN remnants. The corresponding solution for a strong ultra-relativistic blast 
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wave was obtained by Blandford & McKcc (1976). For recent reviews of astrophysical 



shock models see Ostriker & McKcc (1995) and Bisnovatyi-Kogan & Silich (1995) and 
references therein. The basic requirement for the realization of the instant shock solution 
is a short-duration injection of energy E into the shock. However there exist possible 
physical situations of permanent injection of energy into the expanding shock which we 
call the injection shock, when a central source has some time- varying luminosity L = L(t). 
More exactly the instant shock solution is not applicable when the duration of energy 
generation by the central source t s is comparable or exceeds the shock expansion time 
£sh) > tsh' The condition t s > i s h is typical for the early stage of a SN explosion, 
powerful wind from stars and non-steady spherical accretion onto compact objects. The 



other example is the injection shock produced in hidden neutrino sources ( Berezinsky fc 



Dokuc' lacv 2001 ) by successive multiple fireballs after numerous neutron star collisions in 
the dense stellar cluster in a galactic nucleus prior to its collapse into a massive black hole. 
Below we derive the extension of the Sedov- Taylor self-similar spherical shock solution to 
the case of varying in time energy injection by the central source of power form E = At k , 



where A and k are constants. The notations and logistics of "Fluid Dynamics" by Landau 



& Lifshitz (1959, Chapter X, §106) are used in the self-similar solution derivation. 

Let us consider a strong expanding spherical shock in an ideal gas with polytropic 
index (Poisson parameter) 



7 = c p /c v = const, 



(1) 



where c p and c v are the gas heat capacities under constant pressure and volume respec- 
tively. All values on the forward side of the shock discontinuity surface (non-perturbed 
gas side) are designated by index 1, e.g. pi, pi, and behind the discontinuity surface 
(shock cavity side) by index 2, e.g. p2, P2- In a strong shock the pressure behind the 
shock p2 far exceeds the pressure in the non-perturbed gas p\ . The precise definition for 
a strong shock (determined from the shock adiabat) is pz/pi ^> (7 + l)/(7 — 1) and is 
similar to condition u\ >> c s , where u% is the shock expansion velocity with respect 
to the non-moving (non-perturbed) gas and c s = (jp/p) 1 ^ 2 is the sound speed in the 
non-moving gas. The following relations are valid for a strong shock discontinuity: 



l>2 



-Ul, p 2 



7+1 



-Pi, P2 



-piu\. 



7+1 7 — 1 7+I 

These relations are the (outer) boundary conditions for our problem. 



(2) 



2. Self-similar solution ansatz 

The ansatz for self-similar expansion of a strong spherical shock is in the observation that 
gas motion after the shock is determined by only two independent parameters: the initial 
gas density p\ = const and the total shock energy E = const because for a strong shock 
we may put p\ = 0. From these two parameters and from two independent variables, 



V. I. Dokuchaev: Self-similar shock with energy injection 3 
radius r and time t it is possible to construct the only non-dimensional combination 



r (pi/Et 2 ) 1/5 (see e. g. Landau fc Lifshitzj |1959|, |5edov| [1969|, |Stanyukovich| [1969D. As a 



result the gas motion reveals self-similar behaviour when different spherical parts of the 
gas after the shock evolve under a constant value of this non-dimensional combination. 
The law for the shock radius evolution would be 



R = R(t) = (3 ^—j , (3) 

where the constant (3 itself is determined from the exact solution. The corresponding 
velocity of shock expansion with respect to the non-perturbed gas is u\ — dR/dt. The 
gas motion behind the shock would be governed by the non-dimensional self-similar 
variable 

t-m (4) 

which is the relative radius with respect to the instant shock radius R(t). On the surface 
of shock discontinuity £ = 1. 

A helpful hint for finding the generalization of instant shock solution is in the following 
observation: The discussed self-similarity survives if the shock energy E is the power-law 
function of time, E — At k , with A — const and k = const. This is because the self- 
similar variable £ is the power-law function of time. Putting E — At retains the power 
law dependence of the variable £ on time. The only complication is in the changing the 
power index in the definition of £ from 2/5 to (2 + k)/5. Now the case k = corresponds 
to the known instant shock solution and, for example the case k = 1 would correspond 
to the shock with a constant luminosity of the shock source, A = L = const. In the 
following we find the corresponding self-similar solution of fluid equations in a general 
case of arbitrary k. The shock radius in general case evolves as 

, 1/5 

The velocity of shock expansion is 
dR (2 + k)R 

Ul = ^ = — St" ■ (6) 
The non-dimensional self-similar variable is defined as 

C R(t) f3\AJ t(a+*)/B' {l > 
3. Fluid equations 

We use the set of fluid equations for describing the shocked gas motion, as follows. The 
continuity equation in spherical coordinates: 

at or r 



R = R(t)=p[ — \ tV+ k V 5 . (5) 
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The momentum conservation (Euler) equation: 

dv dv 1 dp 
at or p or 

The entropy equation: 



ds 
di 



ds 
dr 



s, 



(9) 



(10) 



where s = s(r,t) is the local entropy source (rate of local entropy production per unit 
mass). The self-similar solution with energy injection (at fc ^ 0) demands the existence 
of some non-hydrodynamic "radiative" mechanism for local energy (and entropy) supply 
from the central source into the shocked gas. In other words it is assumed that additional 
energy is radiatively pumped into the shock from the central source. For the realization 
of self-similar motion the entropy source s(r,t) must be properly tuned. The required 
entropy source term is calculated from the other fluid equations by assuming the self- 
similar motion (Eq. (p5|)). In the particular case of instant shock (k = 0) the entropy 
source term equals zero and the gas motion is adiabatic. 

The local energy equation (energy conservation law) is written as 



0_ 

Of 



2 1 



-div 



pv(w + 



pTs, 



(11) 



where the gas internal energy density e = c v T = of / [7(7 — 1)], w = c p T = 0^/(7 — 1) 
is the gas enthalpy density, T is the gas temperature and c s is the sound velocity. The 
integral form of the energy conservation equation (pd|) is 



0_ 

dt 



dV = - 



pv(w + —) 



di 



pTsdV. 



From the last equation we find the rate of energy injection into the shock 



F- 9 



/ V X 



dV = 



pTsdV, 



(12) 



(13) 



where R = R(t) is the outer shock radius. We assume that this energy injection into the 
shock is supplied by the central source luminosity L = E = kAt k ~ x with the properly 
tuned local entropy source s = s(r, t). 



4. Integral of self-similar motion 

The crucial step for finding the strong spherical shock self-similar solution is the guess of 
the integral of self-similar motion derived from the energy balance equation for a chosen 
sphere expanding in a self-similar manner by the law £ = r/R(t) — const. This law 
of motion determines according to Eqs. (^|) and (m) the corresponding radial expansion 
velocity of this chosen sphere: 
d 



dt 



(14) 



We derive the required integral of self-similar motion in full analogy with Landau fc 



Lifshitz ( 1959| ) procedure for adiabatic (fc = 0) case. The leaking of energy through the 
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sphere surface 4irr 2 of gas with local velocity v and density p during the small time- 
interval At is 

v 2 

pv{w + —)Airr 2 At (15) 

The last equation contains the enthalpy w — e + p/p, where the gas internal energy 
density e = c v T = c 2 /[7(7 — 1)], because the gas produces some work under expansion. 
On the other hand the volume of the sphere increases during this time-interval At by 
the value 47rr 2 Atw n . Inside this volume there is gas with energy 
v 2 

pv n (e+ — )^r 2 At. (16) 

The central source pumps into this volume by a nonadiabatic process the additional 
energy 

pv n TAs4irr 2 At. (17) 

where As = sAt. This additional energy is proportional to the square of At, and At is 
a small time-interval. From the energy balance equation 

2 2 

pv(w + — )4nr 2 At = [pv n (e + —) + pv n TAs]4Trr 2 At (18) 

in the thin spherical shell of radius r < R(t), thickness Ar = v n At and volume 47rr 2 Ar 
in the limit At — > we obtain finally the integral of self-similar motion 

2 2 

v{w+ V -)=v n {e+ V -). (19) 

The left hand side and right hand side of this integral of motion is simply a different 
representation of the same energy flux and does not contain directly the local entropy 
source s — s(r, t). Nevertheless this entropy source is taken into account in Eq. (|h]) indi- 
rectly through connections between the thermodynamic values w and e and the entropy 
relation s — ln(p/p 7 ) for an ideal gas. As a result the required integral of self-similar 
motion ( |l9| ) has the same functional form for both the adiabatic k = an nonadiabatic 
k cases. 

5. Solution of self-similar equations 

We will find the self-similar solution for strong spherical shock with energy injection 
by using three equations: the continuity equation (|§|), momentum equation (^|) and the 
integral of motion ( |l9| ) . The entropy equation ( |l0| ) will be used for the determination of 
the entropy rate source s corresponding to a suitable injection of energy from the central 
source. 

Let us define the non-dimensional variables V — V((,), G — G(£), Z = Z(£) instead 
of gas velocity v, density p and sound velocity c s : 



v = — V; p = piG- c 2 s 



{2 + k)r 
57 



i 2 

Z. (20) 
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According to Eq. (||) on the shock discontinuity surface, i.e. at £ = 1, these non- 
dimensional functions take the values 

V(D = ^; G(1) = X±{; *(1) = ^=£. (2D 
7+1 7-1 (7+1) 



The integral of motion from Eq. ( |19[ ) in non-dimensional variables takes the form 

2-(3 + 7 )y + 2 7 F 2 f , 

v(y-i)(7V-i) 1 ; 

and does not depend on k. With this integral of motion the continuity equation (^) takes 
the form 

x dlnG dV , s 

^iw-dtrr™ (23) 

Correspondingly the Euler equation (^) can be written in the form 

rflnG 2-2(7+2)V+(3 + 27+7 2 )y 2 -2 7 y 3 dV _ 2{5+^[(2 + fc)(y-2)+ 7 (fc-3)]} 
dlnf ( 7 -l)(l-F)l/ 2 ( 7 y-l) dln£~ (7— l)(2+fc)(l— V)V ' ( ^ 

Eqs. ((H), ([||) and ( ^4|) for functions Z, V, and G define completely the motion of the 

shocked gas. By using these three equations and the expression s — c v ln(p/ p 1 ) for the gas 

entropy we now calculate from Eq. ( |Io|) the required entropy source inside the shocked 

gas 

= 7^(7-1)*- (25) 

This entropy source is needed for the implementation of self-similar motion of the shocked 
gas with energy injection in accordance with the law of local energy conservation ( p+| ) or 
the law of energy input into the shock . We will find below from the exact solution 
that 2/(7 + 1) < V < I/7. So ds/dt = in the adiabatic case of k — and ds/dt > in 
the non-adiabatic case of k > 0. See Fig. [|for the radial entropy rate profile for the case 
of k — 1. Meanwhile the entropy source becomes negative, s < 0, at k < 0. It seems that 
a physically reasonable self-similar solution does not exist for the k < case. Resolving 
the system of Eqs. (|4|) and ( p3|) with respect to derivatives dV/dln^ and dlnG/dln£ 
we obtain correspondingly: 

dlnV = (jV - 1){10 - 2[2(2 + k) - 7 (fc - 3)]V + (3 7 - 1)(2 + k)V 2 } 

rflnC (2 + jfe)(l-V)[2- 2(7 + l)V + 7 ( 7 +1)1/ 2 ] 1 j 

and 

rflnG _ y{2(l-y)[l + (2 7 -7)y+47y 2 ]+fc[6-8( 7 +l)y+(7+27+57 2 )F 2 -47y 3 ]} 
dln£ ~ (2 + fc)(l - l/) 2 [2 - 2(7 + T)V + 7(7 + 1)V 2 } ' ( ' 

Now it is possible to integrate the last two equations separately by using the boundary 

conditions from Eq. (^l|). The integration of Eq. ( pif ) is simple but the rather tedious. 

We find the solution for V(£) in parametric form: 



2 10 + ^(2 + fc)(3 7 -l)+2 7 (fc-3)-4(2+fc)] 
17 J 2( 7 -l)(7-7 + 2fc(3 + 7 )] 



a 1 



( 7 +i)( 7 y-i) 



X 



7-1 

(7+l)Vl a2 J (V^D + A)[VD(-y + 1) + 7 2 (fc - 3) + 5 7 (1 +k)- 4(2 + fc)] 1 ^ 



(V# - A)[\+D( 7 + 1) - 7 2 (fc - 3) - 5 7 (1 + k) + 4(2 + fc)] 



(28) 
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in the case of D > and 

£ = I( 7 + 1)2 10 + ^[ y ( 2 + fc )( 3 7 - 1) + 2 7 (fc - 3) - 4(2 + fc)] ^ " 



2(7-l)[7-7 + 2fc(3 + 7)] 



(7+l)V 



( 7 + l)(7^-l) 
7-1 



exp 



2ai 



arctan W 



(29) 



in the case of D < respectively. In these equations the numerical constants 
A = V(3 7 - 1)(2 + fc) + fc( 7 - 2) - 3 7 - 4; 

D = 7 2 (fc - 3) 2 + 2(2 + fc)(9 + 2k) - 2 7 (2 + fc)(9 + 2k); 
(3 7 - 1)(2 + k)[2 - V(7 + 1)]V=S 



W = 
and 

oi = 

a 2 = 
a3 = 

CI4 = 



a 5 {(2 + fc)[(3 7 - 1)V - 2] + 7 (fc - 3)} - £>( 7 + 1) ' 



(7-2)[3(12-7 7 +137 2 ) + (27+l)(3 7 - l)fc 2 ] + fc(33 7 +167 2 + 197 3 -32) 



(30) 



10(27+1) (37-1) 



(2 + fc) 

5 ; 

7-1 



27+1' 

7 [7 + fc + 7(6fc - 13)] -12_--fc 
10(27 + 1)(37 - 1) ' 
a 5 = 4(2 + k) - 7[5 - 3 7 + (5 + 7)*;]. 



(31) 

From Eq. ( |28| ) and ( p9| ) it follows that V(0) = 1/7 is independent of k. 

Similarly we find the solution for G(£) in parametric form by integrating of Eq. (|27j). 
Case D > 0: 



1 



1-V 



(7+l)(7V-l) 



7-1 



(7+1) 1 



10+y[^(2 + fc)(37-l) + 2 7 (fc-3)-4(2+fc)] 
2(7- l)[7-7 + 2fc(3 + 7 )] 



G= 



[ (y/D - A)[VD(j + 1) - 7 2 (fc - 3) - 57(1 + k) + 4(2 + fc)] 
Case D < 0: 



(\AD + A)[V5(7 + 1) + 7 2 (fc - 3) + 57(1 + fc) - 4(2 + fc)] 



x 77 



(32) 



G = 



1 



(7 + IX7V-I) 



1 - V 
(7 + 1) 



7-1 



exp 



2a e 



-D 



arctan W x 



2 10 + F[F(2 + fc)(3 7 - 1) + 2 7 (fc - 3) - 4(2 + fc)] 
2(7 - l)[7-7 + 2fc(3 + 7 )] 



(33) 



where 



« 7 



«8 



3(7 + 3){2 + fc - 7 [2 + 7 (fc - 3) + fc]} 
(3 7 -l)(2 7 +l) 



2 7 + l' 

3(7 2 + 1) 



(2 7 +l)(3 7 -l)' 
6. The calculation of constant (3 from the energy integral 



(34) 



The constant f3 which appears in a self-similar variable £ can be found from the shock 
total energy: 



E 



R / v 2 



P 



7-1 



4irr 2 dr. 



(35) 
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Table 1. Numerical values of constant /3(7, k) for shocks with differ- 
ent 7 and k = (instant shock with E = const), k — 1 (injection 
shock with L = const) and k = (Shock with E = At 2 ). 



7 = 




5/3 


7/5 


4/3 


k = 





1.152 


1.033 


0.994 


k = 


1 


0.929 


0.826 


0.793 


k = 


2 


0.368 


0.288 


0.271 



In non-dimensional variables this can be written as 

Z 



7(7- 1) 



Cd£ = 1. (36) 



By putting into this integral the calculated functions Z, V, G from (22), (|28|) and (|32j) we 
find after numerical integration the required f3 — /3(7, k). This calculation of a numerical 
constant (3 completes the finding of generalized self-similar shock solution. See Table |l] 
for some examples of numerical values of j3 = /3(7, k) 

7. Shock with permanent energy injection 

The particular case of the generalized shock solution for astrophysical application is a 
permanent energy injection shock produced by a continuous pumping of energy into the 
shock from the central source of constant luminosity, L = const. It corresponds to the 
case of k = 1 and E = Lt oc t in contrast to the usually applied Sedov- Taylor shock 
with E = const, which is produced by a instant implosion of energy into the shock. 
The possible astrophysical implications of injection shock solutions with k 7^ 0, that 
is with a central energy source varying in time, are the early phase of SN explosion, 
rarefied bubbles in the interstellar medium after SN explosions, strong wind from stars 
and young pulsars, non-steady spherical outflow from accreting black holes and dense 
stellar clusters near collapse with frequent neutron star collisions. Different analytical 
and numerical approaches were applied by Fallc ( 1975| ), Castor ct al. ( 1975 ), Weaver et 



|aL| ( |1977|) to the modeling of permanent energy injection shocks in the case of stellar 
winds and on the interstellar medium and interstellar bubbles. 

The external radius of the expanding k — 1 injection shock evolves with time ac- 
cording to Eq. (|) as R(t) = f3(A/p 1 ) 1 / 5 t( 2+k '>/ 5 oc i 3 / 5 , where L = A is the luminosity 
of the central source and pi = const is the density of the ambient gas medium. So the 
injection shock expands faster than the instant Sedov- Taylor shock (k=0, R oc t 2 ^ 5 ). This 
is because of a constant pumping of energy into the shock, L — const. The corresponding 
velocity of injection shock expansion is ui(t) = dR/dt = (R/t)(k + 2)/5 oc t~ 2 ^ 5 . From 
the last two equations it is very clear the physical meaning and the difference between 
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the instant shock (k=0) and injection shock (k=l) cases respectively is very clear: 

E = [(5/2) 2 /3(7,0)~ 5 ]i? 3 (pw 2 ) = const; (37) 
L = [(5/3) 3 /3( 7 , l)~ 5 ]i? 2 u(pu 2 ) = const. (38) 

In these expressions R 3 ~ V is shock volume and R 2 ~ S is the shock surface. At 
k = we have the expansion law for instant shock which corresponds to constant energy 
carried by the swept out gas. The expansion law for a k = 1 injection shock corresponds 
to constant energy flux (or constant source luminosity) carried by the swept out gas. 

The discussed strong shock solution is valid only in the region where the shock ex- 
pansion velocity c s <C u(R) <C c, where c s is the sound speed in the ambient gas. The 
expanding strong shock becomes weak and disappears when its expansion velocity drops 
below the the sound speed c s . The maximum radius of the expanding strong shock i? s h 
is obtained from the equality u(-R s h) = c s by using Eqs. (^) and (^|): 

R sh = -=— ■ (39) 

V 5c s J pi 

The maximum time of a strong shock expansion is 



t s \\ — 



2 + k \ 5 A 



P - 



(40) 



5c s ) pi 

The minimal radius of Newtonian motion of the shocked gas, which is called the Sedov 
length Ig, is defined by equality u(Js) — c. The corresponding expression for the Sedov 
length, Is — i? s h(cs/c)( 2+fe ^ 3-fc ) <C R s h, is reproduced from Eqs. ( p9| ) by substituting c 
for c s . So the region of applicability of the strong shock solution is Is <C r <C i? s h- 

8. Conclusion 

The derived solution is the generalization of the Sedov- Taylor self-similar strong spherical 
shock solution for the case of an energy injection from the central source of form E = At k , 
where A and k are constants. The power-law ansatz E = At k only is enough for deriving 
the scaling law for shock radius (^|) and shock expansion velocity (g) accurate to within 
the numerical constant /3(7, k) ~ 1 without knowing the exact solution. The numerical 
value of this constant (see Table |l|) can be calculated from equation Eq. (|3^) only after 
the complete solving of the self-similar problem. The special case of k = corresponds to 
the known Sedov- Taylor solution, while the case k = 1 corresponds to permanent energy 
injection into the shock by a central source of constant luminosity. The cases with k < — 1 
seem to be nonphysical due to the total energy divergence at t — > 0. 

The self-similar hydrodynamic flow in the nonadiabatic k ^ case exists only under 



the self-consistency condition (|25| ) for the local entropy input. In other words the self- 
similar behavior of an expanding shock in the nonadiabatic k ^ case is realised only 
under the appropriate tuning of local entropy (energy) source according to Eq. (|25|). 
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This is the auxiliary physical condition which supposes some radiative mechanism for 
sustained energy supply from the central source into the shocked gas (which depends on 
the detailed properties of the cental engine, radiative transfer, gas composition etc). It 
can be seen from Figs. [I] and ^that the main part of the energy is injected near the outer 
boundary of the shock at £ > 0.8, i. e. at the same place where the shocked gas is mainly 
gathered. The similarities of the profiles for density and entropy rate are in favor of the 
principal realization of the required tuning of the local energy injection mechanism if the 
central source radiation absorption would be proportional to the gas density. 

The self-similar shock solution with energy injection may be applied to the modeling 
of astrophysical objects in which duration of central source activity is longer than shock 
expansion time, such as the early phase of SN explosion, strong wind from stars and 
young pulsars, non-steady spherical outflow from black holes and collapsing dense stellar 
clusters with numerous neutron star collisions. 

Acknowledgements. I am grateful to V. S. Berezinsky and B. I. Hnatyk for useful discussions. 
This work was supported in part by the INTAS through grant 99-1065. 
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Fig. 1. Shock density profile p/p2 = G(£,)(l ~ l)/(7+ 1) depending on the relative radius 
of the shock £ = r/R(t) for 7 = 5/3 and 4/3; k — corresponds to Sedov- Taylor case 
with the shock energy E = const; k = 1 corresponds to the shock with a luminosity of 
the central source L = const. 
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Fig. 2. Shock pressure profile p(r) j-pi = C 2 ^— depending on the relative radius 
of the shock £ = r/R{t) for 7 = 5/3 and 4/3. 
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Fig. 3. Shock velocity profile v/v2 — £V(£)(7 + l)/2 depending on the relative radius of 
the shock £ = r/R{t) for 7 = 5/3 and 4/3. 




Fig. 4. Relative entropy rate (ds / dt) / (c v / 1) profile, where ds/dt from Eq. ( |25| ) 
on the relative radius of the shock £ = r/R(t) for 7 = 4/3 and 5/3, and fc = 



, depending 
0.1 and 1. 



